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Abstract 

The Nordstrom- Vlasov system is a Lorentz invariant model for a self- 
gravitating coUisionless gas. We establish suitable a-priori-bounds on the 
solutions of this system, which together with energy estimates and the 
■ smoothing effect of "momentum averaging" yield the existence of global 

weak solutions to the corresponding initial value problem. In the process 
we improve the continuation criterion for classical solutions which was 
derived recently. The weak solutions are shown to preserve mass. 

Keywords: Nordstrom scalar theory of gravitation; Vlasov equation; global 
weak solutions 
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1 Introduction 

Consider a large ensemble of particles which interact by force fields which they 
create collectively and not by collisions. Such a coUisionless gas is often used as 
a matter model both in plasma physics and in astrophysics. If the particles in 
the gas interact by electromagnetic fields the dynamics of the ensemble is de- 
scribed by the Vlasov-Maxwell system. If the particles interact by gravitational 
forces one obtains the Vlasov-Poisson system in the non-relativistic case or the 
Einstein- Vlasov system in the relativistic case. We refer to 00 for back- 
ground information on these systems. In the present investigation we consider a 
self-gravitating coUisionless gas where gravity is described by Nordstrom's scalar 
theory JOI ■ We write the Nordstrom- Vlasov system in the formulation of ^ : 

dU^^A^^i^, (1.1) 
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^,it,x)^ ffit,x,p)^M=, (1.2) 

J yi+p2 

dtf+P-V.,f- [(^0)p+(l+p2)-i/2v,0] .Vp/ = 4/^0. (1.3) 

Here tsR, x^p^M.^ stand for time, position, and momentum, f = f{t,x,p), (j) = 
(j){t,x)^ and 

Vl+P^ 

S is the free-transport operator, and p denotes the relativistic velocity of a 
particle with momentum p. Units are chosen such that the mass of each par- 
ticle, the gravitational constant, and the speed of light are equal to unity. A 
solution (/,0) of this system is interpreted as follows: The space-time is a four- 
dimensional Lorentzian manifold with a conformally flat metric which, in the 
coordinates {t,x), takes the form 

5^, = e2^diag(-l,l,l,l). 

The particle density on the mass shell in this metric is e~'^'^ f{t,x,e'^p), but it 
is more convenient to work with / and </> as the dynamic variables as long as 
it is kept in mind that / itself is not the particle density. More details on the 
derivation of this system are given in P] where its steady states are investigated. 

Like Nordstrom's theory, the system is not a physically correct model. Nev- 
ertheless, there are good reasons for studying it: The Nordstrom- Vlasov system 
is a Lorentz invariant model for a self-gravitating gas which has the correct 
Newtonian limit, namely the Vlasov-Poisson system, cf. The system is much 
simpler than the physically correct but notoriously difficult Einstein- Vlasov sys- 
tem, and yet it captures some typical relativistic effects, such us the propagation 
of gravitational waves. From a more mathematical point of view, the hope is 
that by investigating this system one can learn more about the whole class of 
non-linear partial differential equations in kinetic theory to which this model 
belongs. For instance, one of the most celebrated results in kinetic theory is 
the existence of global weak solutions for the Vlasov-Maxwell system due to 
R. DiPerna and P.-L. Lions 0. It is the purpose of the present paper to inves- 
tigate this solution concept for the Nordstrom- Vlasov system. 

We explain how the paper proceeds and how global weak solutions to the 
Nordstrom- Vlasov system are obtained. First one needs to establish suitable a- 
priori-bounds for solutions. One set of such bounds derives from conservation of 
energy. In addition, for related systems like Vlasov-Poisson or Vlasov-Maxwell 
the fact that / is constant along characteristics and that the characteristic flow 
preserves measure yields bounds on the L*-norms of f(t), l<q<oo. However, 
these properties do not hold in the present situation. In the next section we 
overcome this diflaculty and derive a-priori-bounds for the L'-norms of f{t). 
This is a necessary prerequisite for global existence of weak solutions, and at 
the same time improves the continuation criterion for classical solutions which 
was derived in ^ . It is interesting to note that the argument would not work if 
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the sign in the field equation were reversed. In the third section we turn to the 
construction of weak solutions. First a suitable regularization of the system is 
introduced such that the regularized system does have global (classical) solutions 
which satisfy the same a-priori-bounds. In |S] the Vlasov-Maxwell system was 
regularized by making the Maxwell equations parabolic. Here we smooth the 
right hand side of the field equation ; an analogous approach was followed 
in [BJ m ini ■ The resulting system is "closer" to the original one than in [5] , in 
particular, it remains time reversible which allows us to avoid various technical 
difficulties in what follows. Along a sequence of solutions to regularized systems 
the a-priori-bounds are shown to hold uniformly. In order to pass to the limit 
in the non-linear term in the Vlasov equation we use the smoothing effect due 
to "velocity averaging" [7|, which we now prefer to call "momentum averaging". 
As opposed to |5] we use this tool only to pass to the limit in the Vlasov equation 
and not for the moments of / such as /x, which are dealt with directly. This 
yields another simplification of the proof. Notice that the Vlasov equation 
can be rewritten in the form 

dtf+p-^.f^Vp- [((50)p+(l+p2)-i/2v^^^ _^fscj, (1.4) 

which is the form needed to apply "momentum averaging" . Once a global weak 
solution is obtained its properties are of interest. For the related systems it 
is not known whether weak solutions are unique or preserve energy, and the 
same is true in the present situation. But we can show that the weak solutions 
preserve mass, which is non-trivial because as opposed to the Vlasov-Maxwell 
system this quantity is not just the integral of f{t), cf. H2.1|) . For this purpose 
we need to exploit the relativistic nature of the system and to assume that the 
initial datum for / is bounded. This assumption is not made in 5 , but it allows 
us in addition to remove further technical difficulties from the proof, such us 
the use of renormalized solutions. 



2 Conservation laws and L^-estimates 

Formally, solutions of the Nordstrom- Vlasov system conserve mass and energy: 

e""* / fdpdx = const, (2.1) 
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f^l+p^dpdx+-J [(9t0)' + (V,0)2]dT = const. (2.2) 

Again formally, these conservation laws can be obtained by integrating their 
local counterparts 

5tp-fV,-j = 0, (2.3) 
ate-fV^-p = 0, (2.4) 
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where 



p{t,x) 




e(t,x) 




Vl+?f{t,x,p)dp+^iM{t,x)f + ^{V,(f>it,x)f, (2.7) 
pf{t,x,p)dp-dtC^{t,x)V,<j){t,x). (2.8) 



Let us now denote by (/,0) G C^{[0,T] x R^) x C^{[0,T] x M^) a classical so- 
lution of the Nordstrom- Vlasov system on the interval [0,T],r>0, with ini- 



tial data /(0) = /'"eCi(K^), 0(0) = 0j)" e C3(R3), 5t<?!)(0) = < e C2(R3), with 



> 0. In the notation above the subscript c indicates that the functions are 
compactly supported while the subscript b indicates that they are bounded to- 
gether with their derivatives up to the indicated order. For such data a unique, 
classical solution exists at least locally in time, cf. Assuming in addition 
that the initial data have finite energy the conservation laws stated above hold. 

Let X{s)=X{s,t,x,p), P{s) = P{s,t,x,p) denote the characteristics of the 
Vlasov equation, i.e., the solutions of the characteristic system 



satisfying X{t)—x, P{t)—p. By the Vlasov equation l|1.3|) the function e ^"^J 
is constant along these curves, and / can be represented as 



In particular, / and hence also /i are non-negative. We write the function (j) as 



where ^hom is the solution of the homogeneous wave equation with initial data 
(/)™, (/)™, and is the solution of with zero initial data. By Duhamel's 
principle, V'f^O. Hence (I2.9|l and (|2.1U|) imply that 



dx 

ds 
dp 

ds 



(50)p-(1+/)-i/2v,0, 



/(i,x,p) = /"XX(0),P(0))exp [40(t,x) -4C(X(0))] . 



(2.9) 



= 0hom+V' 



(2.10) 



f{t,x,p) < f "(X(0),P(0))exp [40hom(i,a:) -4<(^(0))] , 



(2.11) 



and we have proved the following a-priori-bound: 



Proposition 1 For all te [0,T], 



||/(t)||oo<||r||ooexp[4(||0hom(OI|oo + |ICIIoo)]. 
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Notice that under the assumptions on the initial data for <j) made above 
||'/'hom(i)||oo <C(l + with C depending on the data. 

Combining this result with the one in ^ we obtain the following improve- 
ment on the continuation criterion for classical solutions: 

Proposition 2 Initial data as specified above launch a unique classical solution 
to the Cauchy problem for the Nordstrom- Vlasov system on a maximal time 
interval [0,Tmax[, and if 

sup{\p\:{x,p) esupp f{t), 0<i<rmax}<oo (2.12) 

</ien T,„ax = 00, i.e., the solution is global. 

Proof: By 0) Thm. 1] a unique solution exists and is global, provided (j) re- 
mains bounded on [OjTmaxi and the condition 12.12|l holds. But together with 
Proposition^ (|2.12(l implies that the source term /i in the field equation Hl.l|l is 
bounded. Hence a bound on (jj follows from the assumption on the momentum 
support, and the proposition is established. □ 

It is standard that the estimates above result in a-priori-bounds for quantities 
like fi or p: 

Kt,x) ^ f f{t,x,p)—f==+ [ f{t,X,p)—f== 
J\p\<R yjl+p^ J\p\>R yl+P 



<||/(t)iu/ -^.^+R-^ f^/YT^fdp 

J\v\<R V 1 + P J 



\p\<R V-L+F 



<CR'\\f{t)\l 



.+R'^I V^^fdp=c[\\f{t)\\ooj V^Tp^fdp^ 



1/2 



where for the last step we choose 

1/4 



R=(^J V^fdp/wfmc 



Squaring both sides of the estimate for fi, integrating in x, and using Proposi- 
tion ^ and conservation of energy implies the estimate 

11^(0 ||2<Ce2|l'>^-(*)ll~, i£[o^y]. (2.13) 

Here C denotes a positive constant which depends only on ||/'"||oo, ||'/'o'||oo, and 
the energy of the initial data, and which may change from line to line. Similarly, 
using 



\ 3/4 



p{t,x)=e~^l fdp + e-'^ fdp 
\p\<R J\p\>R 



< We- 



^"/"Ilooe3*i?3 + e-'^i?-iy y^l+^fdp^c(^J ^/l+^fdp^ 
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where we have chosen 



1/4 



in the last step. Elevating this estimate to the power 4/3, integrating in x, and 
observing conservation of energy we conclude that 

IIpWII4/3, IIjWII4/3<C; (2.14) 

note that p dominates j. 

Combining (|2.9|l with Liouville's Theorem one can prove the following 
estimates on the distribution function. Since we will assume initial data 
for / in the construction of weak solutions Proposition ^ and the bound on the 
kinetic energy gives us a bound on any L'^-norm of f(t), and hence the following 
estimates are not used in the rest of the paper. But they may be useful if one 
does not wish to consider bounded initial data for /. 

Propositions For all q>l, '^>3/q~4:, andt^[0,T] we have 

\\e^'/'-'^^fm,^\\e^'^'-'^^°nu<\n^ (2.15) 

lle''^/(Ollg<llrll9CXp[7(||0hom(OI|oo + ||C||oo)+|7lll0hom(OI|oo]. (2.16) 

Proof: For any smooth function Q:M^R, Liouville's Theorem and ()2.9|l 
imply 

'"0(/e-'''^)e3'*dpdx = const; (2.17) 



note that the (a;,p)-divergence of the right hand side of the characteristic system 
equals ~3S(j> and so 

det ^^^^ (0, i , x,p) = exp [3 (/.(t, - 3 (X (0, t , a;,p) )] . 

With the choice Q{z)^z'J, l(TT7jl imphes (|TT^ . Moreover 

g[3/g-4]0 J ^ g70g[3/g-4-7]0homg[3/g-4-7]^ j ^ 

Since [3/g — 4 — 7]?/'>0 this yields 

|jg[3/g-4]0^(^)||_^>g-[3/g+4+l7|]||0.o„IU||g70J||^^ 

and (|TT^ follows. □ 
Although we have formulated the results of this section only for going for- 
ward in time they hold equally well towards the past since the system is time 
reversible. 



6 



3 Global weak solutions 

The purpose of this section is to prove global existence of weak solutions to 
the Nordstrom- Vlasov system. We denote by the Banach space of the 

measurable functions g : ^ R for which the norm 



Il5lll,kin = 



J I V^Ti^\g\dpda 



is finite. 

Theorem For any triple (/'", 0q", 0'") such that for some s > 3/2, 

there exists a global weak solution {f,4>) of the Nordstrom- Vlasov system, more 
precisely, 

feL^{R; (M^ ) ) n L^, (M; (M^) ) , e {R;H\R^)), 

with dt(l),V^(j)eL^{R;L^{R^)), e'f' eH^l^iRxR^), f>0 a. e., and the following 
holds: 

(i) {f,(l)) solves ^^j^-^^) in the sense of distributions. 

(a) The mapping 

F:R^L^{R^)xH\R^)xL^(R^), t^ {f{t),(l)(t),dt<t>{t)) 

is weakly continuous with F{Q) — {f™,(j)™,(j/i). Moreover, for any R>0, 
cf>eC{R;L\BR)). 

(Hi) The energy at any time t is bounded by its initial value, the local con- 
servation law 12. holds in the sense of distributions, and the mass is 
conserved: 



p{t)dx = J p{0)dx for a. a. teR. 

Here Bfi = {x G R^ : |x| < R}, and iJ^(M^) denotes the usual Sobolev spaces. 
The proof proceeds in a number of steps: 

Step 1: The regularized system. 

Let < (5„ e C^(M^) be a moUifier satisfying the following conditions: 



5n{x)^5n{-x), J 6n{x)dx^l, supp (5„ C i3i/„, neN. 

The Nordstrom- Vlasov system is regularized by replacing the right hand side 
of (|nj by —fi-kSn~^Sn, where * denotes convolution with respect to x. The 
regularized system is supplied with regularized initial data 
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where 

are such that (f^) is bounded in and as n^oo, 

f'^^r in 4i„ni^ in i^^ /in^^i" in i/^-^. 

The reason for regularizing /i and the data for in this particular way will 
become obvious in the energy estimate below. The regularized initial value 
problem has global smooth solutions. To see this observe that by H2.16|l , | 1 1 1 
is bounded. Hence for any tigN, pi{t)-k6n*Sn is smooth and bounded in L°°, 
together with all its derivatives, the bounds of course depending on n. Hence 
4> and dt4' are bounded on any compact time interval, together with all their 
spatial derivatives. It is then straight forward to see that a standard iteration 
scheme like the one employed in ^ converges on any time interval to a smooth 
solution of the regularized initial value problem. For more details we refer to the 
proof of the analogous result for the relativistic Vlasov-Klein-Gordon system in 
IHl Sect. 3]. Note that the regularized system remains time reversible so that 
solutions really exist on all of M and not only on [0,oo[. 

In what follows, {fm<t>n) denotes the global smooth solution of the regular- 
ized initial value problem. It satisfies the continuity equation 

9t/3„+Vx-j,i=0, 

where pn = e""^" / fndp, j„ = e""*" / pfndp. 
Step 2: Uniform bounds on [fn,4>n)- 

First we observe that the energy of {fm4>n) is not conserved, but one can prove 
that it is bounded. In fact, a direct computation shows that 




Denote by the solution of the wave equation with the given right hand side 
— /in*(5„ and initial data (f>o^n=9n, 4'T.n = ^n- By uniqueness, 4'n = 4>n*Sn- On 
the other hand, the energy of {fm4>n) is constant, and we conclude that 

\\ut)\\iMn+\\\dtKmi + \\\^Mmi 

= ||/«(i)l|l.ki„+^||5t^„(O*'5«||2 + ^l|V:.0„W*5„||^ 
< \\fnmiMn+\\\dtMml + \\\^.~^nit)\\l 

= C- (3.1) 
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constants denoted by C do not depend on n or t. The bound on dt4>n implies 
that \\4>n{t)\\2<CT, for all T>0 and t<^[-T,T], where Ct depends on T but 
not on n. Hence for any T > 0, 

\\4>nmm<CT,t^[-T,T]. (3.2) 

Let 4>hom,n denote the homogeneous part of the field in the modified system. 
By Sobolev estimates for the solution of the homogeneous wave equation and 
the Sobolev embedding theorem it follows that 

I V^noin,n 

(i)||oo<CT, te[-r,r]. (3.3) 

Hence the analogue of (|2.1(J|) and the fact that /j,„*(5„*5„ is non-negative imply 
that 

< e-^-^*) < e'^''""-"^*) < Ct, i G [-T,T]. 

Hence 

l!5te*"(*^||2<||e*"(t)||oo||9t<^4i)ll2<CT, te[-T,T] 
and similarly for the x-derivatives. Hence for any T > and i? > 

\\e^-\\HH\-T.T[.B,)<CT,R, ||e^"(*)||Hi(B«)<CT,fl, te[-T,T]- (3.4) 

clearly the same estimates hold for e^"^". From Proposition H3.3(l . and the 
fact that the sequence {4>^Qn) uniformly bounded it follows that 

e-^^-fn<C, 0<fn,<CT, te[-T,T]. (3.5) 

Hence by Eqn. H2.13(l we conclude from H3.5(l and l|3.1|) that 

\\f^nm2<CT,te[~T,T]. (3.6) 

By (|2.14|l we conclude that 

\\Pn{t)\U/3, \\jnm4/3<C, t GR. (3.7) 

Step 3: The weak limit. 

Here we use the estimates proved in Step 2 to obtain a weakly convergent 
subsequence of {fn,(t>n) whose limit will be a global weak solution. The repeated 
extraction of suitable subsequences is not reflected in our notation. By H3.1f) . 
and (ISi3)-(E!ZI), there exist 

/ e (R; (R6 ) ) n L^, (M; (M^ )),(/) e Lj- (M; ffi (M3 ) ) 

jleL^,{R;L\R^)), p,jeL-(M;L4/3(K3)) 

such that 

/>0 a. e., at./),V,<^Gi:°"(M;L2(R3)), 
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and up a subsequence and for all T> 0, 



/„-/in L2(]-r,T[xM6)^ 
<^„--0in ffi(]-T,T[xR3), 
dtK^dtC^ in L\]~T,T[xm 



Pn,Jn^P,J in L4/3(]-r,r[xM3), 

By a standard diagonal sequence argument we can choose the subsequence and 
its limit independent of T > 0. Next we observe that by (|3.4|l we can choose 
the subsequence such that e"^" converges weakly in H^{]—T,T[xBji) for any 
T>0 and R>0. Since that space is compactly embedded in — T,r[xi?j^) 
it follows that we can choose 0„ and e*^" to converge strongly in L'^, and by the 
Riesz-Fischer Theorem also pointwise a. e. Hence 

e^.^e"^ in Hl^iRxM.^), 

and the same is true for e^"^. We shall prove now that ji — fi, p = p and j—j 
almost everywhere, where p,p and j are related to {/,</)) by (|1.2|) . (|2.5|) . and 
H2.6|l . respectively. As to /i, we have for any x S Cj?°(M x R'^) and any i? > 0, 

/ /'^^ ^ -ti]xdxdt+ [ [ [ f /^^ xdxdt 



^"^^ I I \ I fn^===-pn]xdxdt+ [ [ [ f^^==xdxdt 



n^oo I I \ /|p|<j^ Vl+P^ / J J J\p\>R 

lim / / / fn^==xdxdt+ [ [ [ f^===xdxdt. 

\p\>R Vl+-P J J J\p\>B. vl+-P 



n — >oo 



The last line can be estimated in modulus by 

;^(sup||/„(i)||l,kin+|i/(t)||l,kin) 

where the constant C is independent of R. Since i? > is arbitrary we conclude 
that /i = / f {l+p'^)~^^'^dp a. e., and 

p= f ,f^^=eL^,{R;L^R^)) with pr.^f^ in L\]-T,T[xR^) 

for any T>0. Finally, if we define cr„ = j fndp so that /9„ = e^"^"cr„ then the 
sequence (ct„) is bounded in L^^{R;L'^^^{M.^)), and hence as for p we can show 
that <7n J fdp in i'^^^C — r,r[xM^). On the other hand we have already seen 
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that e*" ^ e"^ strongly in (] - T, T [ x Bfl) . Hence cr„ = e*" p„ ^ e-^p in the sense 
of distributions, and so e'^p— J fdp a. e. Thus we have proved that 

p = e-'^ I fdpeL°°{R;L^/^{R^)) with Pn^pin L'^^^i]-T,T[xW^) 



for any T> 0. It is obvious that the same assertion holds for j. 

Remark 1. It is at this point that we need the -bound on / and hence on 
the initial data: Since the fact that e'^ £ Lf^^{'M.'^) seems optimal we must have p 
in the dual space L^^^ which would not be true if only an L'-bound for / with 
g<oo, say q = 2, were available. 

The pair {4>,p) satisfies the inhomogeneous wave equation p.l|) in the sense 
of distributions. In fact 



= / I {dt<j)n- A(l)„ + Pn*(>n*Sn)xdxdt 

:>nd^x-(f>n^X + fJ-n{Sn-*^Sn*x))dxdt^ j I {(j)dfx~ 4> + t^-x) dx dt 



for any test function x £ ^^^"(M x R'^). The functions p and j satisfy the con- 
tinuity equation (|2.3|l in the sense of distributions. This implies quite easily 
that there exists some constant to which the mass J p{t)dx is equal for almost 
all t, but we want to show that this constant is equal to the initial mass the 
point being that a-priori we have no continuity of the mass as a function of t. 
Conservation of mass in the sense of the theorem is shown in Step 6 below. The 
assertion on the total energy in item (iii) of the theorem is standard. 

Passing to the limit in the linear part of the Vlasov equation (11.311 is not a 
problem. To complete the proof that {f,4>) is a weak solution of (|l.ll) (|1.3|l . we 
need to take care of the non-linear terms in the Vlasov equation. This is done 
in the next step. 

Step 4- Momentum averaging. 

Let i? > and T > 0, ?A G (R^) with supp C Br, and fix ^ G (M) such that 
0<^<1 and C = l on [-T,T]. Define 

where K = (5(/)„)p+ (1 +p2)-i/2 y^^^. Then we obtain 

cf. Eqn. 1)1. 4|) . By the estimates of Step 3 the sequences ((7q) and (g") are 
bounded in L^(R x x Br). Hence, by .5^ the sequence 



fn{-,-,p)'ip{p)dp 

is bounded in H^/^{RxR^), see also 013. Since for all R'>0 and T>0, 
H^/'^{]-T,T[xBr,) is compactly embedded in L'^{]-T,T[xBr,), we conclude 
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that after extracting a subsequence, which, by a diagonal sequence argument, 
can be chosen independent of r> and R' > 0, 



fn{;-,p)^{p)dp^ / f{;-,p)^{p)dp, Strongly in L^]~T,T[xBn,). 



Using this information we can prove that (/, 0) satisfies the Vlasov equation in 
the form ((131). To see this let x(^C'^\]~T,T[xBr,) and ^peC^\BR). Then 
for the non-linear part in (|1.4|) we have 



[Vp-{fnFn) + fnS<pn]x{t,x)^ip)dpdxdt 

-J fn{jp-Vp1pP+{l+P^y^'^ypll)-ll)p)dp-V:^(l)n xdxdt 

f{-p-Vp^j + 'ip)dpdt(l) 



f{p-Vpi^p+ii+p^)-^/^Vp^-^Pp)dp-y^(f> 



xdxdt; 



note that the weight functions in the p integrals of /„ are compactly supported 
smooth functions as required by the momentum averaging argument. Hence the 
p integrals converge strongly in L^{]~T,T[xBji'), and since dt(j)n and V^^n 
converge weakly in the assertion follows. By a standard density argument, 
we conclude that 11.4|l is satisfied in I?'(RxR^). This completes the proof of 
item (?) of the theorem. 

Remark 2. In the application of the momentum averaging argument the L°°- 
bound on / could have been avoided, using the concept of renormalization, cf. 
[B, £ . However, since we already needed this bound above, cf. Remark 1, we 
prefer to avoid this technical complication at this point. 

Step 5: Continuity in t. 

Using the Vlasov equation in the form (|1.4() for the approximating sequence we 
have, for any test function x G C'c 



fn{t)xdpdx= I I f^xdpdx 



fn {s)\p-^xX+{S(pn)is)x 

is^n)is)p+ (1 +pT'/^y.Ms)) ■ Vpx 



dpdxds. 



The convergence of /„ and 0„ is strong enough to pass to the limit in the right 
hand side of this equation, and dropping n in the right hand side we can use the 
resulting expression to define a time dependent distribution f{t) e P'(M^) which 
obviously is continuous in t with respect to the usual topology of P', satisfies 
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dt 



the initial condition /(0) = /'", and by construction coincides for almost all 
t with f(t). By a density argument the I?'-continuity extends to continuity 
with respect to the weak topology of L^(M^). A similar argument works for (j) 
in which case the continuous representative is defined by considering the first 
order formulation of the wave equation Hl.l|l . 

The stronger continuity assertion for (j) follows from the Arzela-Ascoli theo- 
rem: Since dt(j)n is bounded in L°°(R;i^(R'^)) uniformly in n the sequence (0„) 
is equi-continuous as a sequence of L^(M^)-valued functions on K. Moreover, 
for each ieR, (<^n(i)) is bounded in iJ^(R^) which is compactly embedded in 

Step 6: Conservation of mass. 

Since dtPn + divjn = and |j„| £ Pn we have for every R>0 and t > 0, 

[ Pn{t)dx ^ - f Pn{t)dSx+ [ dtPn{t)dx 

'\x\>R+t J\x\ = R+t J\x\>R+t 

Pn(t)dSx- / div jn{t)dx 

x\=R.+t J\x\>B.+t 

{pn{t)+l^-Jn{t))dSx<0 

x\=B.+t 

where v is the outer unit normal of the domain {|a;| > R + t}. The analogous 
argument works for t <0 and the domain {|a;| > R — t}. Hence 

Pn{t)dx< / p„(0)dx, ieR, i?>0, neN. (3.8) 

|>fl+|tl "'|x|>_R 

We claim that for almost all t £ R, 

p{t)dx^ j p{0)dx. (3.9) 

Let e > be arbitrary. Since p{0) is integrable, then we can choose R>0 such 
that 



p{0)dx<e. 

'\x\>R 

By the convergence of the initial data and H3.8|l we conclude that 



Pn{t)dx< / pn{0)dx<e 

x|>J?,+ |t| J\x\>R. 

for all t e M and all sufficiently large n G N. Let A C M be measurable and 
bounded and denote by \{A) its Lebesgue measure. Then 

/ / p(t)dxdt > p{t)dxdt— lim / / pn{t)dxdt 

JaJ j AJ\x\<R.+ \t\ ^^'^JAJ\x\<R+\t\ 



= lim / / pn(t)dx- I Pn{t)dx dt 

^^°°Ja \J J\x\>R+\t\ J 

> X{A)( [ p{0)dx-e), 
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and for sufficiently large 5' > we have by monotone convergence, 



p{t)dxdt< / / p{t)dxdt + \{A)e 

AJ JaJ\x\<S 

= lim / / pn{t)dxdt + X{A)e 
"^°<'JaJ\x\<s 



< XiA) p{0)dx + e 

Since A was an arbitrary, bounded, measurable subset of M this implies that 
there exists a set C M of measure zero such that 

p{0)dx-e< J p{t)dx< J p{0)dx + e, teR\M^. 

Hence holds on M\UfcgNMi/fc. 

Remark 3. The argument above makes use of the relativistic nature of the 
system, i.e., of the finite propagation speed of particles. We do not know if 
conservation of mass in the sense of (|3.9|l holds without this property. In par- 
ticular the above argument would also establish conservation of charge for the 
relativistic Vlasov-Maxwell system, but not for its non-relativistic version in 
which velocity and momentum of the particles are equal. Note that the latter 
system is the one studied in detail in [^. 
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